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Abstract. We consider the nonlinear Dirac equation, also known as the Soler model: 

idtip = -ia ■ VV> + mPi> - f(tp*0tp)fiil!, ip(x,t) e C N , / € C 2 (R), 

where ay, j = 1, . . . , n, and /3 are N x N Hermitian matrices which satisfy cA — ft 2 — In, ctj/3 + /3ay — 0, 
ctjCtk + ctkCtj = 25jfeJ]v. We study the spectral stability of solitary wave solutions 4>{x)e~ t ^ t . More 
precisely, we study the presence of point eigenvalues in the spectra of linearizations at the small amplitude 
solitary waves in the limit uj — s> m. We consider dimensions n = 1, 2, and 3, and prove that if f(s) = 
s k + 0(s k+1 ), k £ N, with m > and either fc > 3, n = 1, or k > 2, n = 2, or k = 1, n = 3, 
then one positive and one negative eigenvalue are present in the spectrum of the linearizations at these 
solitary waves with uj sufficiently close to uj — m. This shows in particular that these solitary waves are 
linearly unstable. The approach is based on applying the Rayleigh-Schrodinger perturbation theory to 
the nonrelativistic limit of the equation. The results are in formal agreement with the Vakhitov-Kolokolov 
stability criterion. 

Resume. Nous considerons l'equation de Dirac non lineaire, aussi connu comme le modele de Soler: 

id t ip = -ia-ViJj + ml3iP-f(ij*f3ij)l34>, i/j(x,t)eC N , 16 1", / G C* 2 (R). 

Nous etudions la stabilite spectrale des ondes solitaires solutions (f>(x)e~ lujt . Plus precisement, nous 
etudions la presence de valeurs propres dans les spectres des linearisations aux ondes solitaires des 
petites amplitudes dans la limite cj — > m. Nous considerons dimensions n = 1, 2 et 3, et montrons 
que si f(s) — s k + 0(s k+1 ), k £ N, avec m > et soit k > 3, n = 1, ou k > 2, n = 2, ou k = 1, n = 3, 
ensuite une valeur propre positive et une negative sont presents dans le spectre des linearisations a ces 
ondes solitaires lorsque uj est suffisamment proche de uj = m. Ceci montre en particulier que ces ondes 
solitaires sont lineairement instable. L'approche est basee sur l'application de theorie de la perturbation 
de Rayleigh-Schrdinger a la limite non relativiste de l'equation. Les resultats sont en accord formel avec 
le critere de stabilite de Vakhitov-Kolokolov. 



1 Introduction 

A natural simplification of the Dirac-Maxwell system [Gro66j is the nonlinear Dirac equation, such as the 
massive Thirring model [Thi58j with vector-vector self-interaction and the Soler model [Sol70j with scalar- 
scalar self-interaction (known in dimension n — 1 as the massive Gross-Neveu model [G N741 |LG75| V There 
is an enormous body of research devoted to the nonlinear Dirac equation, which we can not cover comprehen- 
sively here. The existence of standing waves in the nonlinear Dirac equation was studied in |Sol70| , [CV86J, 
|Mer88| . and |ES95) . Numerical confirmation of spectral stability of solitary waves of small amplitude is 
contained in [Chu07| . An overview of the well-posedness of the nonlinear Dirac equation in ID is contained 
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in [PellO . The asymptotic stability of small amplitude solitary waves in an external potential has been 
studied in |Bou08[ iPSlOj . 

These models with self-interaction of local type have been receiving a lot of attention in particle physics, 
as well as in the theory of Bose- Einstein condensates |HC09| |MJZ + 10 . The question of stability of solitary 



waves is of utmost importance: perturbations ensure that we only ever encounter stable configurations. 
Recent attempts at asymptotic stability [PS IOl [BC11] rely on the fundamental question of spectral stability: 

Consider the Ansatz ip(x,t) = yf> u (po) + p(x, i))e~ lLJt , with 4> ul (x)e~ lujt a solitary wave solution. 
Let dtp = A u p be the linearized equation for p. Does have eigenvalues in the right half-plane? 



In spite of a very clear picture of the spectral stability for nonlinear Schrodinger and Klein-Gordon 
equations |VK73[ ISha83[ [\Vei85, SS85, GSS87] and many attempts at spectral stability in the nonlinear 
Dirac context (let us mention |AS83| IAS86| ISV86|. IBou08l ICKMSlOj ). in the latter case the question of 
spectral stability is still completely open. Numerical results |BC09| show that in the ID Soler model (cubic 
nonlinearity) all solitary waves are spectrally stable. Let us mention the related results |CP06[ [Chu07| . 

Previous results |Comll| show that the Vakhitov-Kolokolov criterion in the case of the nonlinear Dirac 
equation gives a less definite answer about the spectral stability than in the case of the nonlinear Schrodinger 
equation. All we know is that at the value of uj where <9 w <3(cl>) vanishes, with Q{uj) being the charge of the 
solitary wave 4>u,e~ lLJt , two eigenvalues of the linearized equation collide at A = 0, but we do not know where 
these eigenvalues are located when d^Qiu) 7^ 0. This is in contrast to the ground states of the nonlinear 
Schrodinger equation, when the condition d^Q{uj) > is enough to conclude that one positive and one 
negative eigenvalue of the linearized operator move out of A = along the real axis; see |VK73l[GSS87[ [CP03j. 

In the present paper, we show that if Vakhitov-Kolokolov guarantees linear instability for the system ob- 
tained in the nonrelativistic limit, then the same result is true for the corresponding small amplitude solitary 
waves in the original relativistic system. We use our approach to show that the small amplitude solitary 
wave solutions of the nonlinear Dirac equation are linearly unstable if the Vakhitov-Kolokolov condition says 
so. Let us mention that our results apply to the solitary wave solutions which we obtain from the solitary 
waves of the nonlinear Schrodinger equation in the nonrelativistic limit of the nonlinear Dirac. 

According to our results, the spectrum of the linearization at small amplitude solitary waves in the ID 
nonlinear Dirac equation with cubic and quintic nonlinearities has no real eigenvalues; instead, one can prove 
the existence of two purely imaginary eigenvalues. To prove that these solitary waves are indeed spectrally 
stable, one also needs to prove that there are no complex eigenvalues with Re A > 0; we already have partial 
results which we will publish elsewhere. 

Our approach is based on the idea that the family of real eigenvalues of the linearization at a solitary 
wave of the nonlinear Dirac equation bifurcating from A = is a deformed family of eigenvalues of the 
corresponding linearization of the nonlinear Schrodinger equation. The model and the main results are 
described in Section [51 The necessary constructions in the context of the nonlinear Schrodinger equation are 
presented in Section[3] The existence, asymptotics, and linearization of solitary waves of the nonlinear Dirac 
equation is covered in Section |4] The statement of Theorem 12.11 for uj < m follows from Proposition 15. 2\ 
which we prove using the Rayleigh-Schrodinger perturbation theory. 



2 Main result 

We consider the nonlinear Dirac equation 

id t ip=-ia- Vip + m[3iP - f(ijj*/3iP)/3ip, 4>{x,t)eC N , x e R", (2.1) 

with m > and ip* being the Hermitian conjugate of ip. We assume that the nonlinearity f(s) is smooth 
and real- valued, and that 

/(0) = 0. 
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Above, a ■ V = ajj^-, and the Hermitian matrices aj and (3 are chosen so that 



(-ia-V+j8m) =(-A + m 2 )/jv, 
where Jjv is the N x N unit matrix. That is, aj and /3 are to satisfy 

ajOtk + akaj = 25jklN, P 2 = In] Oijfi + (3a,j 



0. 



(2.2) 



The generalized massive Gross- Neveu model, or, in the terminology of [CKMS10J, the scalar-scalar case with 
k E N, corresponds to the nonlinearity f(s) = s k . 

According to the Dirac-Pauli theorem (cf. |Dir281 lvdW32l IPau36| and |Tha921 Lemma 2.25]), the par- 
ticular choice of the matrices aj and /3 does not matter. In terms of the standard Pauli matrices, 



0"! 



r 

1 0) ' 



^2 



-i 

1 



CT3 = 



1 
-II ' 



(2.3) 



which satisfy tJi(jj + <Jj<Ji = 28ijl2, we take: 

n = 1 : = 2, a 



-o- 2j 



n = 2: N = 2,aj = aj, l<j<2, 



n = 3 : = 4, aj ■ = ' ° ° 2 



and in all cases 



N/2 





a 3 







i < i < 3, 



-I 



N/2 



For a large class of nonlinearities /(s), there are solitary wave solutions of the form 
4){x,t) = (f> u {x)e- iut , cj> u £H\R,C N ), \u\<m, 
where in dimensions n = 1, 2, and 3 respectively, 

'1 



(2.4) 



<t>w(x) 



v(x, (jj) 

u(x, lu) 



v{r, uj) 
ie^ufr, lu) 





cos 



v(r, u>) 



(2.5) 



For n = 1, w(x, w) can be taken positive and even, and u(x, uj) real- valued and odd; under these conditions, 
the solitary wave (f>u,(x) is unique (see Section |4~T1 for details). 

For n = 2 and n = 3, v{r, uj) and u(r, uj) are real- valued, radially-symmetric functions, (r, 0) are standard 
polar coordinates in R 2 , and (r,8,(f>) are standard spherical coordinates in R 3 . The existence of solitary 
waves of this form is proved, for example, in [CV86!. The particular solutions we consider here, however, are 
those constructed for uj close to m, as outlined in Section [4.21 from the nonrelativistic limit uj — > m. 

Due to the U(l)-invariance, for solutions to (|2.ip the value of the charge functional 



qwo= / mx,t)\ 2 dx 

is formally conserved. For brevity, we also denote by Q(uj) the charge of the solitary wave 4> ul (x)e~' lujt : 

Q(u) = f |^(s)| 2 daj. (2.6) 



We are interested in the spectral stability of linearization of (|2.1| at a solitary wave solution (|2.4I 
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Theorem 2.1. Let n < 3. Assume that f(s) = s k , where k £ N satisfies k > 3 for n = 1, or k > 2 for 
n = 2, or fc = 1 /or rt = 3. TTien t/ie solitary wave solutions <\> u) e~ lMt to (12. ip described above are linearly 
unstable for we/, where I is the largest interval with sup/ = m such that d^Qicj) does not vanish on I. 
More precisely, let be the linearization of the nonlinear Dirac equation at a solitary wave <p UJ (x)e~ lu)t . 
Then there are eigenvalues 

±A W e a p (A u) ), A w > for lu el, X u = 0(m - u). 

See Figured! 





i(m + lu) 




i(m — ui) 


— A w 


A w 



Figure 1: Main result: The point spectrum of the linearization of the nonlinear Dirac equation in W 1 , n < 3, with 
f(s) = s k , k > 4 — n, at a solitary wave with uj < m contains two nonzero real eigenvalues, ±A W , with A„ = 0(m — cj). 
See Theorem 12.11 Also plotted on this picture is the essential spectrum, with the edges at A = ±i(m — u>) and with 
the embedded threshold points (branch points of the dispersion relation) at A = ±i(m + ui). 

Remark 2.2. Theorem 12. II extends easily to nonlinearities / € C 2 (M.) of the form 

f(s) =as k + 0(s k+1 ), a>0. 

Remark 2.3. We only need to prove the linear instability for u < m. Then, by |Comll| . the positive and 
negative eigenvalues remain trapped on the real axis, not being able to collide at A = for lu 6 / = (loq, to) 
as long as dujQ(uj) does not vanish on / (this nonvanishing is a sufficient condition for the dimension of 
the generalized null space to remain equal to two). These eigenvalues cannot leave into the complex plane, 
either, since they are simple, while the spectrum of the operator is symmetric with respect to the real and 
imaginary axes. 

Remark 2.4. We do not know what happens at to — inf /: as u drops below inf /, it could be that either 
the pair of real eigenvalues, having collided at A = when to — mil, turn into a pair of purely imaginary 
eigenvalues (linear instability disappears), or instead two purely imaginary eigenvalues, having met at A = 0, 
turn into the second pair of real eigenvalues (linear instability persists). 

Remark 2.5. Theorem 12.11 is in formal agreement with the Vakhitov-Kolokolov stability criterion [VK73J, 
since for uj < m one has Q'{ui) > for k > 4 — n. Let us mention that the sign of the instability criterion, 
Q'(ui) > differs from jVK73j because of their writing the solitary waves in the form ip(x)e +lut . 

Remark 2.6. We expect that in the ID case with k — 1 the small amplitude solitary waves are spectrally 
stable; we will prove this elsewhere. 

We also expect that in the ID case with k = 2 ("quintic nonlinearity", f (ip* flip) f3ip = 0(\ip\ 5 )) the small 
solitary wave solutions of the nonlinear Dirac equation in ID are spectrally stable. For the corresponding 
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nonlinear Schrodinger equation (quintic nonlineaxity in ID), the charge is constant, thus the zero eigenvalue 
of a linearized operator is always of higher algebraic multiplicity. For the Dirac equation, this degeneracy 
is "resolved": the charge is now a decaying function for w < to (with nonzero limit as u -> to), suggesting 
that there are two purely imaginary eigenvalues iA^ in the spectrum of A u , with A w = o(to — u>), but no 
eigenvalues with nonzero real part. 

Remark 2.7. Let us notice that in the 3D case for the cubic nonlineaxity f(s) = s (this is the original 
Soler model from |Sol70j), based on the numerical evidence from |Sol70i [AS83| . one expects that the charge 
Q(to) has a local minimum at u = 0.936, suggesting that the solitary waves with 0.936 < u < 1 are linearly 
unstable, but then at u) = 0.936 the real eigenvalues collide at A = 0, and there are no nonzero real eigenvalues 
in the spectrum for uj < 0.936. 

Remark 2.8. We can not rule out the possibility that the eigenvalues with nonzero real part could bifurcate 
directly from the imaginary axis into the complex plane. Such a mechanism is absent for the nonlinear 
Schrodinger equation linearized at a solitary wave, for which the point eigenvalues always remain on the 
real or imaginary axes. At present, though, we do not have examples of such bifurcations in the context of 
nonlinear Dirac equation linearized at a solitary wave. 

Remark 2.9. For n — 3, we only consider the case k = 1, and we do not consider dimensions n > 3. This 
is because of the fact that the equation — Au + u — \u\ p ^ 1 u in R ra has nontrivial solutions in _ff 1 (IR") if and 
only if 1 < p < 2£±| [BL83, Example 1]. This does not allow us to construct solitary wave solutions to the 
nonlinear Dirac equation in R 3 with k > 2 and in R n , n > 3, with any k G N. 

3 Nonlinear Schrodinger and its solitary waves 

We are going to use the fact that the nonlinear Dirac equation in the nonrelativistic limit coincides with the 
nonlinear Schrodinger equation, 

id t i\> = — — A^-|V| 2 V if>(x,t)eC, xeR n , n=l, 2, 3. (3.1) 
2to 

3.1 Solitary waves 

The properties of solitary wave solutions 

^(x,t) = ^(x)e-^ t , cb^eH^W 1 ), 
with amplitude 4>uj(x) therefore satisfying the equation 

A6(x) - \6\ 2k d> = web, lei". (3.2) 

2 TO 

are well-known |Str771lBL83) . For any k £ N when n < 2 and for k = 1 when n = 3, and for any cj E (— oo, 0), 
there is a unique positive, radially symmetric solution 

4> u {x) = MM) > °> 

which decays exponentially. This family of solitary waves, known as the ground states, is generated by 
scalings of a single amplitude function: 

4> u (x) = \u\*SF(y/2m\u\x), uj<0, (3.3) 

where F(x) — F(\x\) > solves 

- AF - F 2k+1 = -F, x e R™. (3.4) 
In one space dimension (n = 1), for k 6 N, there is the explicit formula 

F(x) - ( k + 1 
\ cosh kx 
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3.2 Linearization at a solitary wave 

To derive the linearization of the nonlinear Schrodinger equation (|3.ip at a solitary wave ip(x, t) — (p u (x)e~ ZUJt , 
we use the Ansatz 

and arrive at the linearized equation 



<9 t p = jl(u;)p, p0,*) = 



Bep(x,t) 
hnp(x, t) 



where 



1 

-1 



l_ 

o u 



with l± self-adjoint Schrodinger operators 

1 



l_(w) 



A + |<^| 2fc -o;, l+(w) = l_(w)-2fc|&,| 



2m 



2 A- 



(3.5) 
(3.6) 

(3.7) 



|x|) we take here are radially symmetric, we may consider 



Since the solitary wave amplitudes 4>u(x) = 
the operators 

Vad) l±,rad '■= I, V± restricted to radially symmetric functions. 



(3.8) 



The linear stability theory of NLS ground states is well understood, and can be summarized, in terms of 
their charge (|2.6|) as follows: 



Lemma 3.1 (Vakhitov-Kolokolov stability criterion |VK73| ). For the linearization (|3.5p at a ground state 
solitary wave (f> L0 (x)e~' lult , there are real nonzero eigenvalues ±A G o~d(jl)> A > 0, if and only if ^jQ(^) > 
at this value of u. If so, then ±A G e>d(jl rac ;) are simple eigenvalues, and moreover ker l+ )roe j = {0}. 

Using dS3D, 

Q{uj)= I \<f >LJ {x)\ 2 dx = \uj\ 1/k [ F 2 {y/2m\u\x) dx = CM 1/fe_n/2 , w < 0, (3.9) 



where C = / Rn F 2 (V2~^y) dy > 0. We see from ([SI]) that for cj < 0, Q'(uj) < for fc < 2/ra, Q'(w) = for 
k = 2/n, and Q'(w) > for k > 2/n. Thus: 

Lemma 3.2. Lei n < 3. If k G N satisfies k > A — n, then o~ p ()l r ad) 3 {i^}; / or some A > 0, and m 
particular the NLS ground states are linearly unstable. 



4 Nonlinear Dirac and its solitary waves 

Solitary waves are solutions to (|2.1|) of the form 

V>(z,t) = ^(*)e _<ut ; ^ e ff^R.C^), w £ I 
and as such the amplitude 4>uj(x) must satisfy 

= — ■ Vt/i + mftip — f (tp* pip) /3ip , x G 

4.1 Solitary waves in one dimension 



(4.1) 



We first give a simple demonstration of the existence and uniqueness of solitary waves in one dimension, 
following the article [BC09] , and allowing for more general nonlinearities /(s). 



(j 



Lemma 4.1. Let /(O) = 0. Denote g(s) — m — f(s), and let G(s) be the antiderivative of g(s) such that 
G(0) = 0. Assume that there is ujq < to such that for given uj € (wo, m) there exists T^ > such that 

ujT^ = G(r u ), uj ^ g(r u ), and ujs < G{s) for s € (0,r u ). (4.2) 

Then there is a solitary wave solution ip(x,t) — <j) u (x)e~ tLJt to (12. ip . where 



v(x, Uj) 

u(x, uj) 



(4.3) 



with both v and u real-valued, belonging to 7J 1 (R) as functions of x, v being even and u odd. 
More precisely, for x £ M and uj € (uJo,m), let us define 3^(x,uj) and W(x,uj) by 

JT = v 2 - u 2 , W = vu. (4.4) 

Then 5£{x,uS) is the unique positive symmetric solution to 

d 2 %T = -d x (-2G(3£) 2 + 2uj 2 ,T 2 ), lim &(x,uj)=0 ) (4.5) 

x— >±oc 

and &(x,uj) = — j^d x ^~(x,uj). This solution satisfies &~(0,uj) = r u . 
Proof. From (14.111 . we obtain: 



ujv = d x u + g(\v\ 2 — \u\ 2 )v, 
uju = —d x v — g(\v\ 2 — |u| 2 )u. 



(4.6) 



Assuming that both v and u are real- valued (this will be justified once we have found real- valued v and u), 
we can rewrite (|4.6[) as the following Hamiltonian system: 



d x u = ujv — g(v 2 — u 2 )v = d v h(v, u), 
—d x v = uju + g(v 2 — u 2 )u = d u h(v, u) 



where the Hamiltonian h(v,u) is given by 

h{v,u) = ^{v 2 +u 2 )- l -G{v 2 -u 2 ). (4.8) 

The solitary wave with a particular uj £ (u)q, m) corresponds to a trajectory of this Hamiltonian system such 
that 

lim v(x,uj) = lim u(x,u>) = 0, 

X— >±0O X— f±00 

hence \va\ x ^± 00 S£ = 0. Since G(s) satisfies G(0) = 0, we conclude that h(v(x), u(x)) = 0, which leads to 

uj(v 2 +u 2 ) = G{v 2 -u 2 ). (4.9) 

We conclude from (|4.9[) that solitary waves may only correspond to \uj\ < to, uj ^= 0. 
The functions 2£{x,ui) and W[x,uj) introduced in (|4.4I) are to solve 

dJP = -{v 2 + u 2 )g($r) + ujSC = -±G{3£)g{&) + wX, 1 ' 

and to have the asymptotic behavior lirnui^oo S£(x) = 0, limui^oo &{x) = 0. In the second equation in 
(|4.10j) . we used the relation (|4.9|) . The system (|4. 10|) can be written as the following equation on 3£: 

d 2 x % = -ds-(-2G(,r) 2 + 2uj 2 ^ 2 ) = ±{G{X)g{X) - uj 2 3£). (4.11) 

This equation describes a particle in the potential — 2G(s) 2 + 2ui 2 s 2 . The condition (|4.2[) is needed for the 
existence of the turning point of the zero energy trajectory in this potential at the point s = r u . □ 
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4.2 Solitary waves in the nonrelativistic limit 

In dimensions n = 1, 2 and 3, we consider solitary wave amplitudes <f> u (x) of the forms given in (|2.5[) . 
Substituting these into the nonlinear Dirac equation (|2.1[) , a straightforward calculation results in the system 



ujv = d r u + ^-u + mv- f(v 2 - u 2 )v . , 

um = — d r v — mu + f(v 2 — u 2 )u 

for the pair of real-valued functions v — v(r,Lu), u — u(r,ui). Notice that equation (|4.12p includes the 
n = 1-dimensional case (|4.6p if we interpret r = x. 
Recalling that f(s) = s', we arrive at 

TL — 1 

ho — m)v = d r u H u — fv, (oj + m)u = —d r v + fu, 

r 

with 

/:= (v 2 -u 2 ) k . 

For the nonrelativistic limit, we set 

to 2 — uj 2 — e 2 , < e <C to, 

and rescale v(r, u>) and u(r,u>) as 

v(r,u>) — e*V(er,e), u(r,u)) = e 1+ ~U(er,e). 

Then V, U should satisfy 

(w - m)eiv = e 2 +^(d R U + *j^U) - e^ k fV 
(lu + m)e 1+ iu = -e 1+ id R V + e 1+ i fU 

where R — er denotes the "rescaled variable". Using uj = m — ^e 2 + 0(e 4 ), and taking into account that 

/ = ( e !y2 _ e *+lu 2 ) k = e 2 V 2k + e 4 0(U 2k + V 2k ), 

we re- write the system as 

(-2^7 + °( (2 )) V = d « U + ^ U - v2k+1 + < 2 0{(U 2k + V 2k )\V\) 
(2m + 0(e 2 ))[/ = -d R V + e 2 V 2k U + e 4 0{{U 2k + V 2k )\U\) ' 1 6} 

The rescaled system (|4.13l) has an obvious limit. Formally (for now) setting V(r) = lim e ^o V(r, e )i U(r) = 
lim e ^.o U{r, e), we arrive at 

1 „ n — 1 a 

V = d R U + U-V 2k+1 , 2mU = -d R V. (4.14) 

2m R 

Substituting the second equation into the first yields 



This equation for V(r) is precisely the equation p. 21) for NLS solitary wave amplitudes W with lu = — 
Thus we let V(r) be the (unique) NLS ground state: 

V(r) := (2m)~&F(r), U(r) := -(2m)"^" 1 F'(r), (4.15) 

with F(r) the unique positive spherically symmetric solution to (|3.4[) . We can use this nonrelativistic limit 
to construct nonlinear Dirac solitary waves for e 2 = m 2 — lu 2 <C m 2 : 
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Lemma 4.2. There is u>q < m such that for ui = \J m 2 — e 2 € (u>o,m), there are solutions of (|4.12p of the 
form 



v(r, uj) 



V(er) + V(er) , u(r,uj) = e 1+ ^ U{er) + U(er) , \\V\\ H 2 + \\U\\ H 2 = 0(e 2 ). 



Remark 4.3. In the one-dimensional case n = 1, since the solitary waves are unique (up to symmetries), 
it follows that these asymptotics describe every solitary wave for w close to 1, or equivalently every small 
amplitude solitary wave. 

Proof. The argument parallels that of ( iu;i()K|. where the (more general) nonlinearity f(s) — \s\ 2e , < 9 < 2 
is considered for n = 3. Writing 

V(R, e) = V(R) + V(R, e), V(R, e) = V(R) + V(R, e), 

and subtracting equations (|4.13[) and (|4.14l) . we arrive at 

-±V + 0(e 2 )V 
Zm 



n - 1 , ~ 



= ( Q R + —f^) u ~ ( 2fc + 1)V 2 *V + 0(\V\ 2k+1 V 2 + \V\ 2k+1 ) + e 2 0((U 2k + V 2k )\V\) 
R 

2mU + 0(e 2 )U = -d R V + e 2 V 2k U + e 4 0((U 2k + V 2k )\U\), 
which, setting 

E(R,e) : 

we may re-write as the system 

HS = O m (e 2 ) + 0(e 2 \E\ + \~\ 2 + | ~| 2fc+1 ), 

where 



V(R,e) 
U(R,e) 



Since 



H 

V 



■ -L. + (2k + l)V 2k -(d R + *j±) 
8r 2m 



€ kerH 



£ekerl + , 77 



2m 



Out, 



and kerl +ra( i = {0} (Cf. definitions (|3.7[) . (|3.8p ). we see that kerH = {0}. It then follows from the fact 
that is bounded from L 2 (ir,(C) to H 2 (R n ,C), that H 1 is bounded from H^(R n ,C 2 ) to iJ 2 (R",C 2 ) 

(here L 2 , are the corresponding subspaces of spherically symmetric functions). Hence 



S = H- 1 {O m (e 2 ) + 0(e 2 \E\ + |S| 2 + |~| 2fc+1 )} , 

and since < C||S||#2 (recall n < 3), we arrive easily at 

\\R.H.S.@M\\ m < C\\O m (e 2 ) + 0(e 2 \E\ + \~\ 2 + \~\ 2k+1 ) 
<C{e 2 + e\\E\\ H2 + \\E\\ 2 H2 + \\E\\ 2k ^}, 



(4.16) 



tf 1 



and so we see that for small enough e, the map on the r.h.s. of (14. 16)) maps the ball of radius e in H 2 into 
itself. A similar estimate shows that this map is a contraction, and hence has a unique fixed point H in this 
ball. Finally, we see from (|4~T6f that ||S|| ff 2 = 0(e 2 ). □ 
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ie^* 2 (r,t) 



, i/)(x,t) = 





cos 



(5.1) 



5 Linear instability of small amplitude solitary waves 

Our first observation here is that on spinor fields of the form 

*i(r,t) 
i * 9 < r '*) (e^sin^, 

in dimensions n = 1, n = 2, and n = 3 respectively, the nonlinear Dirac equation (|2.ip reduces to the system 

i&tf i - (9 r + ^)* 2 + m*i - /(|*i| 2 - |* 2 | 2 )*i , , s 

id t ^> 2 = -8r*l ~ + /(|*l| 2 - |* 2 | 2 )*2 1 ' 

(with the convention r = x for n = 1). The solitary waves considered here lie in this class of fields, and so 
to prove their instability it suffices to show that they are unstable as solutions of 



5.1 Linearization at a solitary wave 

To derive the linearization of system (|5.2[) at a solitary wave 
the Ansatz 



where 



" *i(r,t) " 




v( 








. * 2 (r,t) _ 


■( 


u( 


r,uj) _ 


+ 


_ P2(r,t) 




p(r,t) 




Px{r,t) 
_ P2{r,t) _ 


e C 2 . 



we consider solutions in the form of 
~ iut , (5.3) 



Inserting this Ansatz into system (|5.2[) with /(s) = s k and recalling that u(r, uS) and v(r,ui) are real-valued, 
we find that the linearized system for p is 



id t p - 



m - lu — (v 2 — u 2 ) k d r + IL -^- 

—d r —m — L} + (v 2 — u 2 ) k 



1 „.2\fc-l 



p-2k(v z -u z ) 



v —uv 

„,2 



Rep. (5.4) 



We note that the above equation is K- linear but not C- linear, due to the presence of Rc p = 
rewrite equation (|5.4I) in terms of Rep € R 2 and Imp g IR 2 : 



Repi 
Rcp 2 



where J corresponds to 



Rep 
Imp 



J 



= JL(w) 



h 
-h 



Rep 
Imp 



We 



(5.5) 



and the 4x4 matrix operator L(w) is defined by 

L(«) = 

where, writing 
we have 

L_H = 



L_(w) 



/ := /(O0-) = (« 2 (^) - «V,^)) fe , /' := /'(0<M, 



'm-w-/ a r + ^ ' 

— d r —m — u) + / 



m — ui— (v 2 — u 2 ) k d r + IL ^- 

—d r -m - lo + (v 2 — u 2 ) k 



(5.6) 
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and 



L + H = L_H-2/' 



v —uv 
—uv u 2 

2\k 1Uf„,2 



m — uj — (v 2 — u 2 ) k — 2k(v 2 — u 2 ) k 1 v 2 
-d r + 2k{v 2 - u 2 ) h - l uv 



d r + 2=1 + 2k{v 2 - u 



2\fc-l„ 



-m — uj + (v 2 — u 2 ) k — 2k(v 2 — u 2 ) k 1 u 2 



(5.7) 



Let us remind the reader that v — v(r,uj) and u = u{r,u) in (|5.6p ~ (|5.7p both depend on uj. 

Then equation (I5.5[) which describes the linearization of the reduced system (|5.2[) at the solitary wave 
<i>u,e~ %UJt , takes the form 



dt 



Rep 
Im p 



JL(w) 



"Rep" 




L_(w)~ 




Rep 


Im p 




-L+(w) 




Im p 



For the sake of completeness we record here the essential spectrum of the linearized operator: 
Lemma 5.1. cr ess (JL(w)) = iR\i(-l + \w\, 1 - 

Proof. The proof follows from noticing that, due to the exponential spatial decay of v(r,u>), u(r,ui) and the 

Weyl theorem, a ess (JL(w)) = er ess (J(D OT - uj)), with D m = ^ m ® , D m = i<7 2 (3 r + a 3 m. At the same 

time, since -D 2 ^ = —A + m 2 , a ess (D m ) = R\(— m, m), while J commutes with D m and cr(J) = {±i}, one 
concludes that 

a ess (J(D m - uj)) = a ess (i(D m - uj)) U o- ess (—i(D m - uj)) = iR\i(—m + \uj\,m - \uj\). 

□ 



5.2 Unstable eigenvalue of JL for u < m 

Proposition 5.2. Let fc > 3 for n = 1, or k > 2 /or n — 2, or k — 1 /or n = 3. There is loq < m (which 
is on n and k) such that for uj G (wqj to ) i/iere are two families of eigenvalues 



±A W G 0p(JL(w)), A w > 0, A w = 0(m-a;). 



Proof. The relation 



L_ 
-L+ 



= A 



, with <p, 6 C 2 , can be written explicitly as follows: 



0. 



-A 







m — lu — f 


r 


j n—l 




<Pi 







-A 


-d r 


—m 


-J+f 




<P2 


m + LU + f + 2fv 2 




- 2f'vu 


-X 









01 


d r - 2f'vu 


m + ui - 


- f + V'u 2 







-A 




^2 



(5.8) 



We divide the first and the third rows by e 2 = m 2 — ui 2 , the second and the fourth rows by e, and substitute 
R = er, <p 2 = e^2, ^2 = £02, to get 









-A 



rn-^-l g + n-1 



(d R + V ) " — 



-m+u+f+2f'v 

d R - -2f'vu m + uj- f + 2f'u 

A 



',,2 



_ A_ 




-R 1 R 

m — uj 





<P1 












02. 



(5.9) 



Anticipating the e — > limit, formally set A = lim -j-, and introduce the matrices 

e-vO e 



k-A = 



-A 

_^ + (2fc + l)V^(i?) + 

2m 

11 



^-V 2fe (i?) 9« + 2=i 
— Or —2m 
-A 







(5.10) 



Ki = diag[l, 0,1,0], K 2 = diag[0, 1,0,1], 
where V(R), the NLS ground state, was introduced in f)4. 15[> . We write (|5.9p in the form 



fx 

$2 
(>1 
02 



eC 4 



(5.11) 



(5.12) 



where W(R, e) is given by 



-J fFi + (2k + l)V 2k + ™->-f-V'" a 






2f'vu 



X 
2m 



2f'vu 



m-uj + / + 2/'u 2 
is a zero order differential operator with L°° coefficients. 
Lemma 5.3. \\W(-, £)||l-(r+,c 4 ^c 4 ) < 0{e 2 ). 
Proof. By Lemma 14.21 (and the Sobolev inequality) , one has 

v(r,uj) =ei(V(R) + V(R)), u{r,u) = e 1+ i (U(R) + U(R)), 

Then 

f(v 2 -u 2 ) = e 2 (V 2 -U 2 ) k , 
f(v 2 -u 2 ) = ke 2 -i(V 2 -e 2 U 2 ) k -\ 



y2k _|_ ui — m+f 









-m + u) — f 





(5.13) 



||V"||i« + ||#||i« =0(e 2 )- 



',,2 



1 m — cj — f — 2f'v 
2m e 2 



= 0{e z ) - {V' z - e 2 U 2 ) k - 2kV z (V' z - e z U' z ) K ~ l =-{! + 2k)V~ ZK + L ~> {e z ), 



2/T/2 2 TT 2\k-X 



r2k 



f'vu 



ke 2 UV(V 2 - eU 2 )^ 1 = L ~(e 2 ) 



m - to + / = 0{e 2 ) + e 2 (V 2 - e z U 2 ) k = O l ~ (e z ) 



and 



/V = ke*XJ z (V z - ef/ 2 ) ft = 0(e 4 ) L 
1 ui — m + / 



= +0{e z ) + {V 2 - e 2 U 2 ) k = V 2k + L <*>(e z ), 



2k 



2m 



and the Lemma follows directly from this list of estimates. 
Lemma 5.4. dimker Aa = dimker(jl rQ d — A), where 



□ 



l+,rad 
^ — ,rad 



and where, we recall, 



1 



,rad 



l+.raJ = -^(Sji + ^JSji + ^ + (2* + l)V(R)' k . 

Zm ti, Zm 

Moreover, if n = 1 and k > 3, or n = 2 and k > 2, or if n — 3 and k — 1, there is A > such that 
±A £ crdQH) fl^e simple eigenvalues. 
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Proof. An easy computation shows that 



$ = 



$1 

$2 
4>4 



£ ker Aa 



2ra$ 2 = -clfi*!, 2m$ 4 = -<9 fl 4> 3 , OVod - A) 



$ 3 



= 0, 



and the first statement of the Lemma follows from this observation. The second statement then follows from 
Lemma 13.21 □ 



So we may assume that there is A > such that ±A £ ^(Aa), with eigenfunctions 



ker A ± a 3 $±/ 



±$i 
-J-d R $ 3 

2m u ° 



l+,rad^l = -A$ 3 , l_, rad $ 3 = A$i. 



We will use the Rayleigh-Schrodinger perturbation theory to show that there are A € <7<z( JL) with A f=s ±e 2 A. 
Writing 

H-o.4- " 



A A = JL -AK l! L = 



L ,_ 



with 



±--(2k + l)V™ d R + ^ 
—8r —2m 



self-adjoint, we see that 



A A = L J AKi = FA A F, F := 



1 T/2fc a i 71—1 

^ - V or + — 
—dn —2m 



h 
h 



Hence ker A A is spanned by <i> A := F<I> A . Let Pa denote the orthogonal projection onto $ 
Seeking 7/ and A in the form 

r? = $A + C, C^-$A, A = e 2 (A + /j), 

then (|5.12[) becomes 

AaC = mKi($a + C) + e 2 (A + /i)K 2 ($ A + C) + W(<5> A + (). 
Applying Pa and 1 — Pa to (|5.14p . one has: 

= ,i($ A , Kx(<i> A + C)> + e 2 (A + /i)($ A , K 2 ($ A + 0) + (* A , W(* + 
AaC = (1 - Pa) (a*Ki + e 2 (A + M )K 2 + wj (* A + 0. 



A- 



Next, we claim 
To see this, note that 



($a,Ki$a> = (F$a,Ki$ a ) = 2Re($ 3 ,*i) 



(5.14) 

(5.15) 
(5.16) 

(5.17) 



Now the fact that <I> A 6 ker A A means in particular that L_$3 = A$i. Hence, since L_ is self-adjoint, 

M3 ($ 3 ,£-*3> = A($ 3 ,$l), 

and so Re($ 3 ,$i) = only if (<£> 3 , L_$ 3 ) = 0. As is well-known, since L^V = and V(r) > 0, we have 
L_ > 0. Thus Re(4> 3 , $i) = only if i_$ 3 = 0. This, in turn, would imply either A = or $ A = 0, both 
of which are false. 



13 



Denote by L^(K™,C 4 ) C L 2 (R™,C 4 ) the subspace of spherically symmetric functions. Now using (|5.17p . 

Ia 1 : Range (1 - P A ) -> 



and the existence of the bounded inverse A A 4 : Range (1 — Pa) — > the equations (|5.15p . (|5. 16[) can be 



written as 

fi = M(fi,(), C = Z(fi,0, 
with functions M : R x L 2 r ->• K, Z : M x 4 given by 

M ^< C ) = ~ 7a= I Jh \ N*A, KxC) + e 2 (A + m)K 2 ($ a + C) + W)l , (5.18) 

Z( M , C) = AX X (1 - Pa) (mKi + e 2 (A + M )K 2 + ($ A + 0- (5.19) 
Pick r > 1 such that 

r>2||AX 1 (l-PA)K 1 $ A ||L2. (5.20) 
Lemma 5.5. Consider KxL^ endowed with the metric 

HfoOllr = r|A*| + ||ciU». 

There is u)\ G (ujQ,m) such that for uj G (wi,m) the map 

MxZ: I x -> 1 x L 2 ., (ftO^(%0-%0). ( 5 - 21 ) 

restricted onto the set 

6 £ = {(/i,()elx L 2 (R, C 4 ); IK/i, C)||r < e} C R x L 2 
is an endomorphism and a contraction with respect to \\ ■ ||r- 

Proof. Assuming + ||C||l 2 < e < 1, and using Lemma T5.3[ we have the estimates 

\M(jm, 0\<c (|a*|||C||w + e 2 (i + ImI)(i + IKIUO + \\W\\ L -) < Ce 2 , 
Olli- < ir| Ai |^c^(| M |||ciU-H-e 2 (i^| M |)(i^||cili-)^l|w|| i ^(i-i-||cl| jE .-)) < C7e2 

which show that for all sufficiently small e, M x Z maps B c into itself. Similar estimates show that (M xZ)\ B 
is a contraction in the metric || ■ ||r- d 

According to Lemma 15 . 51 by the contraction mapping theorem, for uj G (u}\ , m) , the map (|5 . 2 1|) has a 
unique fixed point (/xo(w), Co(<^)) e 8 e C i x L 2 . Thus, we have ±e 2 (A + /i (u;)) G er p (JL(u;)), w G (u>i,m), 
with r|/io(w)| < e, finishing the proof of the proposition. □ 



By Remark [231 Proposition 15.21 finishes the proof of Theorem 12. II 

References 

[AS83] A. Alvarez and M. Soler, Phys. Rev. Letters, 50 (1983), pp. 1230-1233. 
[AS86] A. Alvarez and M. Soler, Phys. Rev. D, 34 (1986), pp. 644-645. 

[BC09] G. Berkolaiko and A. Comech, On spectral stability of solitary waves of nonlinear Dirac equation 
on a line, ArXiv e-prints, 0910.0917 (2009). 

[BC11] N. Boussaid and S. Cuccagna, On stability of standing waves of nonlinear Dirac equations, ArXiv 
e-prints, 1103.4452 (2011). 

[BL83] H. Berestycki and P.-L. Lions, Nonlinear scalar field equations. I. Existence of a ground state, Arch. 
Rational Mech. Anal., 82 (1983), pp. 313-345. 



14 



[Bou08] N. Boussaid, On the asymptotic stability of small nonlinear Dirac standing waves in a resonant 
case, SIAM J. Math. Anal., 40 (2008), pp. 1621-1670. 

[Chu07] M. Chugunova, Spectral stability of nonlinear waves in dynamical systems (Doctoral Thesis), Mc- 
Master University, Hamilton, Ontario, Canada, 2007. 

[CKMS10] F. Cooper, A. Khare, B. Mihaila, and A. Saxena, Solitary waves in the nonlinear Dirac equation 
with arbitrary nonlinearity , Phys. Rev. E, 82 (2010), p. 036604. 

[Comll] A. Comech, On the meaning of the Vakhitov-Kolokolov stability criterion for the nonlinear Dirac 
equation, ArXiv e-prints, 1107.1763 (2011). 

[CP03] A. Comech and D. Pelinovsky, Purely nonlinear instability of standing waves with minimal energy, 
Comm. Pure Appl. Math., 56 (2003), pp. 1565-1607. 

[CP06] M. Chugunova and D. Pelinovsky, Block- diagonalization of the symmetric first-order coupled-mode 
system, SIAM J. Appl. Dyn. Syst., 5 (2006), pp. 66-83. 

[CV86] T. Cazenave and L. Vazquez, Existence of localized solutions for a classical nonlinear Dirac field, 
Comm. Math. Phys., 105 (1986), pp. 35-47. 

[Dir28] P. A. M. Dirac, The quantum theory of the electron, Proc. Roy. Soc. A., 117 (1928), pp. 616-624. 

[ES95] M. J. Esteban and fi. Sere, Stationary states of the nonlinear Dirac equation: a variational ap- 
proach, Comm. Math. Phys., 171 (1995), pp. 323-350. 

[GN74] D. J. Gross and A. Neveu, Dynamical symmetry breaking in asymptotically free field theories, Phys. 
Rev. D, 10 (1974), pp. 3235-3253. 

[Gro66] L. Gross, The Cauchy problem for the coupled Maxwell and Dirac equations, Comm. Pure Appl. 
Math., 19 (1966), pp. 1-15. 

[GSS87] M. Grillakis, J. Shatah, and W. Strauss, Stability theory of solitary waves in the presence of sym- 
metry. I, J. Funct. Anal., 74 (1987), pp. 160-197. 

[Gua08] M. Guan, Solitary wave solutions for the nonlinear Dirac equations, ArXiv e-prints, 0812.2273 
(2008). 

[HC09] L. H. Haddad and L. D. Carr, The nonlinear Dirac equation in Bose-Einstein condensates: foun- 
dation and symmetries, Phys. D, 238 (2009), pp. 1413-1421. 

[LG75] S. Y. Lee and A. Gavrielides, Quantization of the localized solutions in two-dimensional field theories 
of massive fermions, Phys. Rev. D, 12 (1975), pp. 3880-3886. 

[Mer88] F. Merle, Existence of stationary states for nonlinear Dirac equations, J. Differential Equations, 74 
(1988), pp. 50-68. 

[MJZ+10] M. Merkl, A. Jacob, F. E. Zimmer, P. Ohberg, and L. Santos, Chiral confinement in quasirela- 
tivistic Bose-Einstein condensates, Phys. Rev. Lett., 104 (2010), p. 073603. 

[Pau36] W. Pauli, Contributions mathematiques a la theorie des matrices de Dirac, Ann. Inst. H. Poincare, 
6 (1936), pp. 109-136. 

[PellO] D. Pelinovsky, Survey on global existence in the nonlinear Dirac equations in one dimension, ArXiv 
e-prints, 1011.5925 (2010). 

[PS10] D. E. Pelinovsky and A. Stefanov, Asymptotic stability of small gap solitons in the nonlinear Dirac 
equations, ArXiv e-prints, 1008.4514 (2010). 

[Sha83] J. Shatah, Stable standing waves of nonlinear Klein-Gordon equations, Comm. Math. Phys., 91 
(1983), pp. 313-327. 



15 



[Sol70] M. Soler, Classical, stable, nonlinear spinor field with positive rest energy, Phys. Rev. D, 1 (1970), 
pp. 2766-2769. 

[SS85] J. Shatah and W. Strauss, Instability of nonlinear bound states, Comm. Math. Phys., 100 (1985), 
pp. 173-190. 

[Str77] W. A. Strauss, Existence of solitary waves in higher dimensions, Comm. Math. Phys., 55 (1977), 
pp. 149-162. 

[SV86] W. A. Strauss and L. Vazquez, Stability under dilations of nonlinear spinor fields, Phys. Rev. D 
(3), 34 (1986), pp. 641-643. 

[Tha92] B. Thaller, The Dirac equation, Texts and Monographs in Physics, Springer- Verlag, Berlin, 1992. 

[Thi58] W. E. Thirring, A soluble relativistic field theory, Ann. Physics, 3 (1958), pp. 91-112. 

[vdW32] B. L. van der Waerden, Die gruppentheoretische Methode in der Quantenmechanik , Springer- Verlag, 
Berlin, 1932. 

[VK73] N. G. Vakhitov and A. A. Kolokolov, Stationary solutions of the wave equation in the medium with 
nonlinearity saturation, Radiophys. Quantum Electron., 16 (1973), pp. 783-789. 

[Wei85] M. I. Weinstein, Modulational stability of ground states of nonlinear Schrddinger equations, SIAM 
J. Math. Anal., 16 (1985), pp. 472-491. 



16 



